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Abstract. The inclusive distributions of gluons and pions for high-energy NN 
collisions are calculated. The results for several unintegrated gluon distribu- 
tions (UGD's) from the literature are compared. We find huge differences in 
Oh. both rapidity and pt of gluons and tt's in NN collisions for different models of 

UGD's. The Karzeev-Levin UGD gives good description of momentum distri- 
bution of charged hadrons at midrapidities. We find, however, that the gluonic 
mechanism discussed does not describe the inclusive distributions of charged 
k^ , particles in the fragmentation region. Some of the missing mechanisms are 

^^ ' calculated with the help of unintegrated parton distributions being solutions 

5^ , of the Kwiecihski CCFM equations. We predict an asymmetry in 7r+ and 7r~ 

production which could be studied by the BRAHMS collaboration at RHIC. 

Keywords: pion production, nucleon-nuclcon collisions, unintegrated parton 
distributions, fragmentation functions, inclusive distributions 
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1. Introduction 

The recent results from RHIC (see e.g. [ITl) have attracted renewed interest in better 
understanding the dynamics of particle production, not only in nuclear collisions. 
Quite different approaches have been used to describe the particle spectra from 
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the nuclear coUisions [|2- The model in Ref.[|ni with an educated guess for UGD 
describes surprisingly well the whole charged particle rapidity distribution by means 
of gluonic mechanisms only. Such a gluonic mechanism would lead to the identical 
production of positively and negatively charged hadrons. The recent results of the 
BRAHMS experiment concerning heavy ion collisions [ W put into question the 
successful description of Ref. [^ . In the light of this experiment, it becomes obvious 
that the large rapidity regions have more complicated flavour structure. 

We discuss the relation between UGD's in hadrons and the inclusive momentum 
distribution of particles produced in hadronic collisions. The results obtained with 
different UGD's [|31 EJ [71 HJ 1^1 are shown and compared. 

In addition to the gg —^ g mechanism we include new mechanisms qfg — > qf and 
gqf — > qf and similar for antiquaries. We present first results based on unintegrated 
parton (gluon, quark, antiquark) distributions obtained by solving a set of coupled 
equations developed recently by Kwiecihski and collaborators. 

2. Inclusive gluon production 

At sufficiently high energy the cross section for inclusive gluon production in hi + 
^2 ^ g can be written in terms of the UGD's "in" both colliding hadrons [ I10| 



dyd^pt iV2 - 1 p2 

Above J-'i and J-2 are UGD's in hadron hi and ft.2, respectively. The longitudinal 
momentum fractions are fixed by kinematics: a;i/2 = -y= • exp(±7;). The argument 
of the running coupling constant is taken as fl^ — max(K^, K^^Pt)- 

Here we shall not discuss the distributions of "produced" gluons, which can be 
found in [^^- Instead we shall discuss what are typical values of xi and X2 in the 
jet (particle) production. Average value < xi > and < X2 >, shown in Fig^ only 
weakly depend on the model of UGD. For y - at the RHIC energy W = 200 GeV 
one tests UGD's at Xg = 10"'^ - 10"^. When \y\ grows one tests more and more 
asymmetric (in xi and X2) configurations. For large \y\ either xi is extremely small 
{xi < 10"'') and 2:2 ^ 1 or xi ^ 1 and X2 is extremely small {x2 < lO""*). These 
are regions of gluon momentum fraction where the UGD's is rather poorly known. 
The approximation used in obtaining UGD's are valid certainly only for x < 0.1. 
^ In order to extrapolate the gluon distribution to Xg —^ I we multiply the gluon 
distributions from the previous section by a factor (1 — Xg)"", where n = 5-7. 

3. From gluon to particle distributions 

In Ref.[l2j it was assumed, based on the concept of local parton-hadron duality, 
that the rapidity distribution of particles is identical to the rapidity distribution of 

^an interesting discussion concerning the kinematics of the standard 2— »2 and the 2—^1 mech- 
anism discussed here can be found in [ 118) 
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Fig. 1. <xi> and <X2> for pt > 0.5 GeV and at W = 200 GeV. 



gluons. In the present approach we follow a different approach which makes use 
of phenomenological fragmentation functions (FF's). For our present exploratory 
study it seems sufficient to assume Of^ = Og. This is equivalent to i^h — r/g — Vg, 
where rjh and rjg are hadron and gluon pseudorapitity, respectively. Then 



j/g = arsinh ( — ^ sinhy,i 



Pt,h 



(2) 



where the transverse mass rrit , 



' pf fi- In order to introduce phenomeno- 



logical FF's one has to define a new kinematical variable. In accord with e~^e~ and 
ep collisions we define a quantity z by the equation Eh — zEg. This leads to the 
relation 

(3) 



Pt,g = J(mt.h,yh) 

z 



where J(jnt^h,yh) is given in Ref. [irij. Now we can write the single particle distri- 
bution in terms of the gluon distribution as follows 



dygCppt^g / dz Dg^hiz,fiD) 



da{-qh,ptji) 
dijhd'^Ptji 
./ ^ x2 f - ^Pt-A da{yg,ptj 

In the present calculation we use only LO FF's from [E] or [ I18j . 



(4) 
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Fig. 2. Charged-pion pseudrapidity distribution at W = 200 GeV for different 
models of UGD's. In this calculation pt^t > 0.2 GeV. The experimental data of the 
UA5 collaboration are taken from [I14|. 



Let us analyze now how the results for pseudorapidity distributions depend on 
the choice of the UGD. In Fig[21l compare pseudorapidity distribution of charged 
pions for different models of UGD's. In this calculation FF from [ 12 has been 
used. 

In contrast to Ref. [El, where the whole pseudorapidity distribution, includ- 
ing fragmentation regions, has been well described in an approach similar to the 
one presented here, in the present approach pions produced from the fragmenta- 
tion of gluons in the gg — > g mechanism populate only midrapidity region, leaving 
room for other mechanisms in the fragmentation regions. These mechanisms in- 
volve quark/ antiquark degrees of freedom or leading protons among others. This 
strongly suggests that the agreement of the result of the gg —^ g approach with the 
PHOBOS distributions [|2| in Ref.[|3] in the true fragmentation region is rather due 
to approximations made in [T than due to correctness of the reaction mechanism. 
In principle, this can be verified experimentally at RHIC by measuring the 7r"'"/7r~ 
ratio in proton-proton scattering as a function of (pseudo)rapidity in possibly broad 
range. The BRAHMS experiment can do it even with the existing apparatus. 

In Figl^l we compare the theoretical transverse momentum distributions of 
charged pions obtained with different gluon distributions with the UAl collabo- 
ration data [El- The best agreement is obtained with the Karzeev-Levin gluon 
distribution. The distribution with the GBW model is much too steep in compar- 
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ison to experimental data. This is probably due to neglecting QCD evolution in [ 
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Fig. 3. Transverse momentum distributions of charged pions at W = 200 GeV for 
BKK1995 FF and different UGD's. The experimental data are taken from [ I15j . 



4. Unintegrated parton distributions from the solution of the 
Kwiecihski CCFM equations 

Many unintegrated gluon distributions in the literature are ad hoc parametrizations 
of different sets of experimental data rather than derived from QCD. Recently Jan 
Kwiecinski and collaborators [ El ^| have shown how to solve so-called CCFM 
equations by introducing unintegrated parton distributions in the space conjugated 
to the transverse momenta [ I16| 



/,(.,»^.=).±/exp(..S) /,(,*, 



fi^) <Ph . 



(5) 



The approach proposed is very convenient to introduce the nonperturbative effects 
like internal (nonperturbative) transverse momentum distributions of partons in nu- 
cleons. It seems reasonable to include the nonperturbative effects in the factorizable 
way 

fCCFMf^ U ,,2^ _ T^np(u\ (a\ 



U{xA^^')^!^'^'''\xA^^')■F^^{h). 
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Fig. 4. Leading-order diagrams for parton production 

In the foUowing, for siniphcity, we use a flavour and a:;-independent form factor 

52 



F"P(6) = F"''(6)=exp 



which describes the nonperturbative effects. 
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(7) 



5. From unintegrated parton distributions to meson produc- 
tion 

The og — > g mechanism considered in the hterature is not the only one possible. In 
Fig01we show two other leading order diagrams. They are potentially important in 
the so-called fragmentation region. The formulae for inclusive quark/antiquark dis- 
tributions are similar to formula (Q and will be given explicitly elsewhere [ I19| . The 
inclusive distributions of hadrons (pions, kaons, etc.) are obtained through a con- 
volution of inclusive distributions of partons and flavour-dependent fragmentation 
functions 



Dg^h{z,tJ?u) 



da(nh,pt,h) 

dr]hCPpt,h 

dcr^g^g(yg,Pt,g) 

dygdPpt.g 



^ 2 x'^'^?/s^'J/fe'J/'Pt,9/) 



/=-3 



dyqjd?pt,q 



dz^- 






Pt,q = ->Pt,h/' 
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Fig. 5. Pseudorapidity distribution of charged pions at W = 200 GeV calculated 
with the CCFM parton distributions. The experimental data are taken from [I15|. 
The thin solid line is the gluon-gluon contribution while the dashed lines repre- 
sent the gluon-(anti)quark and (anti)quark-gluon contributions. In this calculation 
gluon, (anti)quark fragmentation functions from [E] have been used. 
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In FiglSlwe show the distribution in pseudorapidity of charged pions calculated with 
the help of the CCFM parton distributions [[T7| and the Gaussian form factor {TJ 
with bo = 0.5 GeV~^, adjusted to roughly describe the UA5 collaboration data. Now 
both gluon-gluon and (anti)quark-gluon and gluon- (anti)quark fussion processes can 
be included in one consistent framework. As anticipated the missing up to now 
terms are more important in the fragmentation region, although its contribution in 
the central rapidity region is not negligible. More details concerning the calculation 
will be presented elsewhere [119). 

For completeness in Fig|Slwe show transverse momentum distribution of positive 
and negative pions for different incident energies. The presence of diagrams Bi and 
B2 leads to an asymmetry in tt^ and 7r~ production. The higher the incident energy 
the smaller the asymmetry. This is caused by the dominance of diagram A at high 
energies. 
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Fig. 6. Transverse momentum distribution of tt^ (solid) and tt 
in proton-proton collisions for different center-of-mass energies. 

6. Conclusions 



(dashed) mesons 



We have calculated the inclusive distributions of gluons and associated charged 
tt's in the NN collisions through the <?5 — > <? mechanism in the fct-factorization 
approach. The results for several UGD's proposed recently have been compared. 
The results, especially pt,h distributions, obtained with different models of UGD's 
differ considerably. 

Contrary to a recent claim in Ref.[in], we find that the gluonic mechanism 
discussed in the literature does not describe the inclusive spectra of charged particles 
in the fragmentation region, i.e. in the region of large (pseudo)rapidities for any 
UGD from the literature. The existing UGD's lead to the contributions which 
almost exhaust the strength at midrapidities but leave room for other mechanisms 
in the fragmentation regions. We conclude that at presently available energies the 
gluonic mechanism is not the only contribution. 

We propose new mechanisms, neglected so far in the literature, which involve 
also quark/antiquark degrees of freedom and are based on (anti)quark-gluon and 
gluon-(anti)quark fusion processes followed by the subsequent fragmentation. These 
missing mechanisms have been estimated in the approach based on unintagrated 
parton (gluon, quark, antiquark) distributions originating from the solution of a 
set of couple CCFM equations proposed recently by Kwiecihski and coworkers. 
They lead to an asymmetry in the production of the 7r+ and 7r~ mesons. While at 
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lower energies of the order of 20 - 50 GeV such an asymmetry has been observed 
experimentally, it was not yet studied at RHIC. The BRAHMS collaboration at 
RHIC has a potential to study such asymmetries. 
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